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We investigate the variation with light quark mass of the mass of the nucleon as well as the masses 
of the mesons commonly used in a one-boson-exchange model of the nucleon-nucleon force. Care is 
taken to evaluate the meson mass shifts at the kinematic point relevant to that problem. Using these 
results, the corresponding changes in the energy of the 1 So anti-bound state, the binding energies 
of the deuteron, triton and selected finite nuclei are evaluated using a one-boson exchange model. 
The results are discussed in the context of possible corrections to the standard scenario for big bang 
nucleosynthesis in the case where, as suggested by recent observations of quasar absorption spectra, 
the quark masses may have changed over the age of the Universe. 



(N 

>: 

ov 

O' 
O' 
cn; 

00' 

o: 

(N 



I. INTRODUCTION 

In the last decade there has been considerable interest 
in the possibility that the fundamental "constants" of 
Nature may actually change with time Although it 
remains controversial, there is growing evidence that the 
fine structure constant may have varied by an amount 
of order a few parts in 10 -5 over a period of 5-10 billion 
years [2 5] . It has even been suggested that this variation 
may have a dipole structure as we look back in different 
directions @. Although this possible variation is quite 
small, within the framework of most attempts at grand 
unification, a variation of a implies considerably larger 
percentage changes in quantities such as Aqcd and in 
the quark masses pHSl- For example, in Ref. pj it was 
shown that the variation 5m q /m q would be of order 38 
times that of 5a/ a. 

In the light of these developments it is very natural to 
ask what other signatures there may be for such changes. 
These may, for example, be the consequent changes in 
hadron masses or magnetic moments |10i4l2l |. Indeed, 
in some cases the level of precision possible in modern 
atomic, molecular and optical physics means that it may 
even be feasible to detect the minute variations expected 
under the hypothesis of linear variation until the present 
day over a period as short as a year [T^-llBl • 

Another consequence of a variation in the parameters 
relevant to hadron structure is the possibility of observ- 
able consequences in big bang nucleosynthesis (BBN) 
or other nuclear phenomena such as the composition 
of the ash of long extinct natural nuclear reactors fisl — 
[Lsl ]. In this context, the effect of quark mass changes on 
the nucleon-nucleon force has been studied in effective 



field theory [19j, |20ll. m ost recently including constraints 
from lattice QCD [5l|, [22[. For the moment these lat- 
tice studies are at too high a quark mass to provide an 
accurate constraint [22T [23| . In an alternative approach 
based on a more traditional model of the nucleon-nucleon 
force, the latest work of Flambaum and Wiringa [24| on 
this topic involved the study of the variation of nuclear 
binding with quark mass using the Argonne potential 
and Schwinger-Dyson estimates of the variation of me- 
son masses. 

In this work we employ a one boson-exchange (OBE) 
model of the nuclear force to calculate the variation with 
changes in the quark mass of the binding energies of se- 
lected finite nuclei as well as the energy of the 1 So anti- 
bound state and the binding energies of the deuteron 
and triton. Apart from its intrinsic interest, this ap- 
proach complements the work of Ref. [24[ and a com- 
parison of the two provides one way to gauge the pos- 
sible model dependence of the variations reported. The 
method used here involves a detailed study of the varia- 
tion of the mass of each of the mesons usually employed 
in a one-boson-exchange (OBE) picture of the nucleon- 
nucleon (NN) force. Care is taken to estimate this shift 
at the relevant kinematic point, not just at the real, on- 
shell meson mass or its pole position. These changes are 
then introduced into the quark-meson coupling (QMC) 
model for some light nuclei and a typical OBE model for 
the two nucleon systems and a Faddeev calculation of the 
triton. 

In section II we examine, in turn, each of the mesons 
do, or, u), p, ir and r\. Section III presents results for fi- 
nite nuclei, while the two nucleon system and triton are 
discussed in section IV. The final section is reserved for 
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some concluding remarks. ™^ 
II. MESON MASSES 



In order to find the variations in the mass of the ex- 
changed bosons in the nuclear interaction (cr, p, and u>) 
with respect to changes in the mass of the pion m v , 
we use three ideas. First we introduce a description of 
the bare mass of the a (mi ) in terms of m^ using the 
Nambu- Jona-Laisino (N JL) model [25|, [26[ . Second, we 
introduce the contribution from the self energies for a, 
p, and uj. Finally, we include these self energies in two 
different ways: for the a the loop diagram is fixed in such 
a way that the total propagator contains a pole on the 
second sheet of the complex e nerg y plane at the posi- 
tion found by Leutwyler et al. [27| : for the p and lo we 
use the chiral fit to partially quenched data from lattice 
QCD developed by Armour et al. [Uj]- Finally, through 
the Gell-Mann-Oakes-Renner (GMOR) relation we relate 
those changes to variations in the quark masses. 



A. VARIATION IN m a WITH m, 

In this work we choose to parametrize the intermedi- 
ate range attraction in the nucleon-nucleon (NN) force 
in terms the exchange of a cr meson, following the tra- 
ditional one-boson-exchange (OBE) approach. Earlier 
work on the effect of changes in quark masses by Flam- 
baum and Wiringa [24j , used explicit two-pion exchange 
for this purpose. Almost certainly the reality is some- 
where in between these extremes and a comparison be- 
tween our results and those of Ref. [24| should serve to 
pin down the uncertainties in this sector of the calcula- 
tion. 

The existence of the a meson has been somewhat con- 
troversial, largely because its width is comparable with 
its mass. However, a careful dispersion relation treat- 
ment using the Roy equation has served to accurately lo- 
cate a pole which can be unambiguously identified with 
the a meson. Of course, because of the large imaginary 
part of the energy of this pole, one cannot easily relate 
the position of the pole to the position of a bump in the 
7T7T cross section. When it comes to the mass of the vir- 
tual a meson exchanged in a OBE NN potential, it is 
a third value that is of interest. Indeed, the invariant 
mass of a meson exchanged in a typical NN interaction 
is very near zero and so we actually need the sigma mass 
for p 2 ~ 0. This is most readily found within an effective 
Lagrangian approach. 

Our model for the a meson involves a "bare" a me- 
son coupled to two pions. When required, the variation 
of the mass of this bare state with quark mass will be 
calculated within the NJL model. In this approach, the 
propagator of the dressed a is described as a bare scalar 
propagator plus an infinite series of contributions of the 
form shown in Fig. [TJ Calling this self-energy the 



FIG. 1. Self-energy contributions for the a meson. 

total propagator can be written as: 
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The pole in A a is the mass of the a resonance. This pole 
was calculated by Leutwyler et al. using the method of 
Roy equations, which is model independent [28|. They 
obtained a pole located in the complex second sheet for 
p at 
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The real part of the position of the resonance, m a — 
441 MeV, is in the range (400 - 1200) MeV given by the 
Particle Data Group (PDG) [H, while its width, T = 544 
MeV, is within the range (600 - 1000) MeV, also from 
the PDG. 

Having a reliable value for the a pole we can find a 
relation that lets us fix such that 



(m^ 0) ) - (ml) ~ 441 - i272 MeV. (3) 

With derivative coupling of the bare a to two pions (con- 
sistent with chiral symmetry), the expression for the itit 
self-energy is found to be: 
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where k represents the pion loop momentum, p the a 
momentum and 7 the aim coupling (initially we took the 
value 70 from Harada, Sannino, and Schechter [29]). We 
are considering all these particles as elementary, so this is 
just an effective theory, and like any other effective theory 
it has to be regularized. The regularization scheme we 
choose is to impose a dipole cut-off (at each vertex) on 
the loop momentum with mass A: 



A 2 



(5) 
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TABLE I. Parameters fixed to reproduce the position of the 
a meson pole (70 = 6.416 x 1CT 3 (MeV^ 1 )). (Am, is the 
deviation of the fitted from the empirical value.) 

7 (x 7o ) A (Mel/) mS 0) (MeV) Am, (MeV) x 10" 5 



4.56 
4.60 
4.70 
4.83 
5.02 
5.27 
5.61 
6.07 



320.00 
330.00 
340.00 
350.00 
360.00 
370.00 
380.00 
390.00 



563.54 
600.23 
639.85 
683.46 
732.00 
790.18 
859.15 
945.43 



1.2 - 1.0? 

1.1 - 0.5i 

1.2 - O.Oi 
1.2 - O.Oi 
0.0 + 0.0i 
0.9 - 1.2? 
0.4 - 0.9i 
0.9-0.9i 



TABLE II. Calculations 
the fractional change of 
the OBEP model to the 



for the coefficient v a , which relates 
the mass of the a meson relevant to 
fractional change in the quark mass. 
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563.54 


-0.145 


2.677 


2.822 


0.089 


600.23 


-0.164 


2.632 


2.796 


0.078 


639.85 


-0.189 


2.576 


2.765 


0.068 


683.46 


-0.220 


2.546 


2.766 


0.060 


732.00 


-0.261 


2.502 


2.763 


0.052 


790.18 


-0.314 


2.451 


2.765 


0.045 


859.15 


-0.389 


2.401 


2.790 


0.038 


945.43 


-0.495 


2.344 


2.839 


0.032 



which is sufficient to ensure convergence. This dipole 
regulator contains simple poles in k (after writing it in 
the form of derivatives with respect to A), which permits 
us to use contour integration over the time component. 
For the remaining integral over the three-momentum 

we rotated k in the complex plane k e 10 , with — ^ < 

9 < 0, to ensure that the imaginary part is located in the 
complex second Riemann sheet. We also performed a 
numerical integration with the help of the routine Nlnte- 
gral of Mathematica. The final value of EJ T (p 2 ) depends 
on two parameters: the regularization mass A, and the 
coupling constant 70. We choose a range of values for A 
such that, after fixing 70 and to reproduce the pole 
position (Eq. ©), m£- 0) varies from 560 to 950 MeV. The 
results are summarized in Table HI where Am a represents 
the deviation of our result for the pole position from that 
of Leutwyler and collaborators. 

We then define m 2 (OBE) = (m { °A - ££ w (0), be- 
cause in an OBEP model for the NN interaction the ex- 
changed boson has nearly zero momentum, m 2 (OBE) 
is a real value because £J„. (0) is real. Thus any variation 
on m a {OBE) with respect to is given by variations 

^ (0) and TZ 
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and using the Gell-Mann-Oakes-Renner (GMOR) rela- 
tion (13: 



Sm a (OBE) 
m a (OBE) 
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We change m^ near the physical value and find the varia- 
which is almost constant for a small change 



tion 



as; (o) 



in m 2 , so we only need to find the slope of the plot 

(0) versus m 2 . We also need the variation of 
with m^ near 140 MeV. For this purpose we used the NJL 
model, which is known to respect the chiral behaviour of 
QCD, including the GMOR relation. The results for all 
the cases in Table U are contained in Table [TTJ 

From Table HH we notice that as A increases (A in- 
creases when m§ grows) the value of £^„. (0) also grows 
and its contribution of S f^i to Sm "g^f E ' > increases. 
However, in practice the variation of the a bare mass 
is much larger and, in total, the larger mi the smaller 
the coefficient v a . What is more, calculations within the 
Quark Meson Coupling (QMC) model tend to favour val- 
ues for m a (OBE) near 550 MeV (see sect. HV[) . 



B. VARIATIONS IN m p AND 
RESPECT TO m„ 



WITH 



In the case of the p meson we have a good deal of 
data taken from lattice calculations in partially quenched 
QCD from the CP-PACS collaboration. Armour et 
al. [32| used this data in an analysis that included the 
leading and next-to-leading non-analytic chiral correc- 
tions to the self-energy to make an extrapolation of the 
mass ra p to the chiral limit (to t w 0). At the physi- 
cal value of m-Tr they found excellent agreement with the 
physical value. 

The relevant self-energy diagrams for the p are given 
in Fig [2] 
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(a) (b) 
FIG. 2. Self-energy contributions for the p-meson. 
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These yield the following expressions: 



yp — _ 

7T7T — 
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6^" 



*r (fc) («* (fc) 



(9) 



Using for m p (OBE), henceforth simply written as m p , the 
value 770 MeV, which is usually used in OBE models, we 
find: 



^ = 0.034^i 



(18) 



yp - 

7TLU — 



12^ 



(k)dk 



o 



- 6.028 and g up 
use 

Uttw (k) 



0.016 MeV" 



where /, 

ularization functions used in the analysis are: 

A 4 



(fc) 



(A 2 + k 2 y 



A 2 + ^-/£ 
(A 2 + fc 2 ) 2 



(10) 
The reg- 

(11) 
(12) 



and we use the approximations: <C m WjP and m p sa 
m u . 

The fit to the partially quenched lattice QCD data for 
the p meson involved a fit of the form: 



nip = \J (a + a2ml + a^m%f + Etot , 



(13) 



where Etot = + E£ w , and the coefficients, <2i, are: 
a = 832.00 MeV, a 2 = 4.94x 10~ 4 MeV -1 , a 4 = -6.10x 
lO" 11 MeV -3 and A = 655.00 MeV (up to errors). At 
the physical pion mass (in full QCD) this yields a value 
of: 



Trip n 778 MeV . 



(14) 



which shows remarkable agreement with the physical 
value, with a shift of only: 



iP hys - 3.7 MeV . 



(15) 



As in the case of the a meson, we consider a one boson 
exchange potential with almost zero momentum transfer, 
so that Hp ~ in the propagator of © (not in the regu- 
lator, because the mass that appears there is the physical 
mass): 



The analysis for the ui meson is closely related to that 
of the p meson. However, the diagrams that contribute 
to the self-energy terms differ because there is no two- 
pion contribution because of G-parity. In addition, 
is 3 x , because there are three possible p — ir charge 
combinations. For the analytic terms in the expansion 
we use the same coefficients (at ) as in the case of the 
p, because the mass difference between them is only of 
order 10 MeV. Then the variation of with respect to 



to 2 near the physical value gives 
which leads to the relation: 



8m u 
A in 2 



0.00096 MeV" 



^ = 0.024^ 



(19) 



where we used the physical mass for the w, m u = 782 
MeV - again because that is the value typically used in a 
OBE potential. (The value obtained at zero momentum 
transfer would be 765 MeV.) 



SUMMARY OF MESON MASS VARIATION 



TABLE III. Coefficients Vi summarising the rate of variation 
of the masses of the mesons used in an OBE description of 
the NN force with respect to quark mass - see Eq. (|20|) . 



Meson 


v (MeV) 


7T 


0.5 


V 


0.012 


cro 


0.089 


<T1 


0.072 


P 


0.034 


Ul 


0.024 



yp — 

7T7T — 



1 



2 



k 4 ul„ (k) dk 

W 3 (fc) 



(16) 



This of course changes the value of m p , now denoted 
77i p (OBE), at the physical pion mass. Indeed, in this 
case it is near 762 MeV. The relation between m p and m 2 
near the physical value is almost linear and it presents 
a slope of Sm " s { °? E) = 0.00135 MeV^ 1 . Following the 
analysis for m a we can relate this change with m q in the 
following way: 



5to p (OBE) 



6m p (OBE)\ 8m, 
J m q ' 



m p (OBE) 5ml 



For the 77, like the pion, we use the GMOR relation to 
calculate the variation within respect to u and d mass. 
In the case of the iso- vector scalar meson, cr±, which has 
negative G-parity and therefore does not couple to two 
pions, we use the NJL model - corresponding to the third 
column and second row of Table UH and Eq. (J5J without 
the self-energy part. For convenience, in Table IIIII we 
summarise the values of za, defined as 



Srrii 5m„ 

= Vi ^ 

TUi m q 

which will be used below. 



(20) 
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III. NUCLEON MASS 

In order to compute the variation of nuclear binding 
energies with quark mass, we also need to know how the 
nucleon mass changes. 

The variation with light quark mass is directly given 
by the so-called ttN sigma commutator 



ovjv = m q (N \qq\ N) = m q 



SlTlN 

5m a 



(21) 



where qq = uu + del. 

The last equality, which gives the information we need, 
follows from the Feynman-Hellmann theorem. A number 
of methods have been used to extract ct„n from pion- 
nucleon scattering data using dispersion relations, but 
the resulting value is still controversial. 

Instead, the most reliable method seems to be to use 
fits to lattice QCD data for mjv as a function of m q [33| . 
These fits, which build the constraints of chiral effective 
field theory, appear to yield very reliable values. We 
take the result of the latest analysis of PACS-CS data by 
Shanahan et. al. [Hj], namely <i t jv = 45 ± 6 MeV. Thus 
we use: 



Stun 
m N 



0.048^. 

m a 



(22) 



IV. 7 LI, 12 C AND 16 O NUCLEI 

To study the effect of the quark mass variation on the 
single-particle energies of 7 Li, 12 C and 16 nuclei, it is 
highly desirable to use a nuclear model based on the 
quark degrees of freedom. The quark-meson coupling 
(QMC) model, which originated with Guichon [37j as a 
description of nuclear matter and was extended and im- 
proved to describe the properties of finite nuclei [H, H^] , 
is ideal for this purpose. The successful features of the 
QMC model applied to various nuclear phenomena and 
hadronic properties in a nuclear medium, are reviewed 
extensively in Ref. |40| . The model has been updated 
to study the properties of hypernuclei [4l| , and neutron 
star structure [42, 4j| , where the quark structure of the 
nucleons and hyperons should play an important role at 
such high density. We calculate the change in the single- 
particle energies of these nuclei versus the current quark 
mass (m q ) and the mass of the nucleon (mjv) using the 
theory presented in Ref. [39| and the meson and nucleon 
mass changes calculated above. 

In Ref. [39] the standard values used to repro- 
duce the nuclear matter saturation properties are, 
(m N , m <r ,m w ,m p )=(939,550,783,770) MeV, with the 
current quark mass m q — 5.0 MeV. For the calculation 
of finite nuclei, the ratio for the a — N coupling constant 
and the mass, (g^ /m a ), was kept constant and fitted to 
the rms charge radius of 40 Ca, r c /,( 40 Ca) = 3.48 fm, by 
adjusting ra a — > fh a = 418 MeV. (Note that the variation 



of m„ at fixed (g^ /m a ) has no effect on the nuclear mat- 
ter properties.) To account for this, we calculate the shift 
in m a from 550 MeV for a given variation of the quark 
mass. This changes the ratio (g^ /m a ) and from that 
new value we deduce the corresponding shift in m CT from 
418 MeV, to be used in the finite nucleus calculation. 
First, with the nucleon mass fixed at mjv = 939 MeV 
and the variations of the meson masses, 8m a ^^ p , evalu- 
ated for quark mass variations of 5m q = ±0.05 and ±0.1 
MeV, we calculate the single-energies in 7 Li, 12 C and 
16 0. Note that the very small differences for the a and uj 
meson mass values used to extract the relation in terms 
of 5m q in III Al and III Bl were neglected. In addition, we 
also calculate the energy per nucleon (E/nucleon). The 
results are given in Table IIVI 

From Table IIVI we see that the absolute values of the 
single-particle binding energies of each nucleus decrease 
as the quark mass increases. This is because an increase 
of the quark mass leads to a significant increase of the 
mass of the a meson and this reduces the attraction aris- 
ing from a meson exchange by more than the repulsion 
associated with the oj decreases. It is interesting to point 
out that a small variation of the quark mass of 0.05 MeV 
is reflected in a change in the single-particle energies of 
order of 0.1 MeV. That is, the impact is appreciable. Fur- 
thermore, we note that the binding energy per nucleon 
for each nucleus decreases linearly as the quark mass in- 
creases. 

Next, we calculate the variation of the single-particle 
energies as the mass of the nucleon is varied. The results 
are given in Table [V] for the same nuclei as in Table IIVI 
As the value of the nucleon mass increases the absolute 
values of the single-particle binding energies also increase. 
This seems to be natural, since the kinetic energy is sup- 
pressed. 

It may be helpful to consider the binding energy per 
nucleon as a function of the quark mass. Based on the 
results given in Tables [TVl and fVl and in (|22p ; we get the 
following relations for each nucleus: 



S\Er u 


| / nucleon 


|Er u | 


/nucleon 


(5|Ei2 C 


| /nucleon 


|El2 C | 


/nucleon 


S |El6Q 


| /nucleon 


|Eie | 


/nucleon 



-2.571^, 
m q 

. Sm q 
m q ' 

. 8m q 
m„ 



= -1.438- 



= -1.082- 



(23) 
(24) 
(25) 



The contributions to the previous coefficients from the 
variation of the exchanged mesons masses are found from 
Table IIVI and from Table |V] we obtain the contribu- 
tion from the nucleon mass. These calculations are sum- 
marised in the following equation: 

S | E, | /nucleon 5m q 

i-r^ i '/ ; y^mesons ± ^Nucleon) ) 

|Ej| /nucleon m q 
with i representing each of the three nuclei we are con- 
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TABLE IV. Single-particle energies (in MeV) for 7 Li, 12 C and 16 nuclei versus quark mass m q (in MeV) calculated in the 
quark-meson coupling (QMC) model [3||. E/nucleon stands for energy per nucleon. The standard value for the quark mass 
use in the QMC model is rn q = 5.00 MeV. 





States m q 


4.90 


4.95 


5.00 


5.05 


5.10 


7 Li 


P 


lSl/2 


-19.0216 


-18.8657 


-18.7089 


-18.5522 


-18.3953 




lf3/2 


-3.5945 


-3.5104 


-3.4267 


-3.3432 


-3.2602 


11 


lSl/2 


-18.3503 


-18.2146 


-18.0777 


-17.9408 


-17.8039 




1^3/2 


-3.2385 


-3.1709 


-3.1036 


-3.0366 


-2.9699 




E/nucleon 


-1.710 


-1.662 


-1.614 


-1.566 


-1.519 


12 c 


P 


lSl/2 


-26.2579 


-26.1110 


-25.9643 


-25.8179 


-25.6716 




1^3/2 


-10.0448 


-9.9445 


-9.8444 


-9.7447 


-9.6453 


11 


lSl/2 


-29.4077 


-29.2558 


-29.1040 


-28.9525 


-28.8011 




1^3/2 


-12.9435 


-12.8379 


-12.7327 


-12.6277 


-12.5231 




E/nucleon 


-4.174 


-4.107 


-4.040 


-3.974 


-3.908 


16 o 


P 


lSl/2 


-29.0713 


-28.9260 


-28.7810 


-28.6362 


-28.4917 




1P3/2 


-13.8605 


-13.7503 


-13.6405 


-13.5309 


-13.4217 




lPl/2 


-12.0635 


-11.9647 


-11.8661 


-11.7679 


-11.6700 


11 


lSl/2 


-33.0861 


-32.9358 


-32.7857 


-32.6358 


-32.4862 




1P3/2 


-17.6266 


-17.5112 


-17.3962 


-17.2815 


-17.1671 




lPl/2 


-15.8151 


-15.7110 


-15.6073 


-15.5038 


-15.4006 




E/nucleon 


-6.109 


-6.035 


-5.961 


-5.888 


-5.815 



sidering, v mesons being described by 

6 | Ej | /nucleon m q 



'mesons 



Sm a 



|Ej| /nucleon' 



and VNudeon by 



_ 8 |Ej| /nucleon m N 

^Nucleon — J ' m i , ; ' 0.U48. 

omjv | Ej | /nucleon 



V. VARIATION IN THE ENERGIES OF THE 
TWO- AND THREE-NUCLEON SYSTEM WITH 
VARIATION IN THE MESON AND NUCLEON 
MASSES 

To examine the variation in the binding energy of the 
deuteron and triton with changes in the meson and nu- 
cleon masses, we need to consider a purely One Bo- 
son Exchange (OBE) model for the nucleon-nucleon in- 
teraction. We choose to employ the OBE potential of 
Bryan-Scott (BS) [44{, which includes the exchange of 
(tt, 77, oo, (Ji, p, a;)-mesons. To avoid the singular nature 
of this potential, BS introduced a monopole regulariza- 
tion scheme that insured that the potential is finite at 
the origin. With a cutoff mass of 1500 MeV this regular- 
ization is shorter in range than the range of the heaviest 



of the bosons included in the potential. As a result, the 
medium range interaction is dominated by the 00 , <j\ fol- 
lowed by the p and u> exchanges. 

Because of the nonlocal nature of the potential (term 
proportional to V 2 ), we have used the method of mo- 
ments [45[ to solve the Schrddinger equation for the bind- 
ing energy of the deuteron and the 1 So amplitude. This 
entails expanding the radial wave function ipe{f) for a 
given angular momentum I as a linear combination of 
Yamaguchi [46| wave functions ipf(r; Pi) with different 
range parameters Pi, i.e. (47j 



Mr) = b t4 Y \ r ->& 



(26) 



where we have taken n = 12 and the Pi are multiples 
of the pion mass. The present choice for the variational 
wave function ensures that the correct long-range behav- 
ior of il>t{r) is that defined by the asymptotic behavior of 

tp£ (r). This in turn is determined by the binding energy 
of the deuteron or the : So anti-bound state. This proce- 
dure reduces the two-body Schrodinger equation to a set 
of 2n homogenous algebraic equations that give us the 
binding energy and the wave function for the deuteron 
to a very good approximation [45L . 

For the ^o, the pole in the scattering amplitude is on 
the second energy sheet, and the analytic continuation of 
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TABLE V. Single-particle energies (in MeV) for 7 Li, 12 C and 16 nuclei versus nucleon mass mjv (in MeV), calculated in the 
quark-meson coupling (QMC) model [3|| . E/nucleon stands for energy per nucleon. The standard value for the nucleon mass 
used in the QMC model is m N = 939.0 MeV. 





States mjv 


938.0 


938.5 


939.0 


939.5 


940.0 


7 Li 


P 


lSl/2 


-18.6768 


-18.6929 


-18.7089 


-18.7249 


-18.7409 




lf>3/2 


-3.4010 


-3.4139 


-3.4267 


-3.4396 


-3.4524 


u 


lSl/2 


-18.0517 


-18.0647 


-18.0777 


-18.0907 


-18.1036 




1P3/2 


-3.0830 


-3.0933 


-3.1036 


-3.1138 


-3.1241 




E/nucleon 


-1.600 


-1.607 


-1.614 


-1.621 


-1.628 


12 c 


P 


lSl/2 


-25.9445 


-25.9544 


-25.9643 


-25.9742 


-25.9841 




1P3/2 


-9.8211 


-9.8328 


-9.8444 


-9.8561 


-9.8677 


11 


lSl/2 


-29.0828 


-29.0934 


-29.1040 


-29.1147 


-29.1253 




1P3/2 


-12.7076 


-12.7202 


-12.7327 


-12.7452 


-12.7577 




E/nucleon 


-4.022 


-4.031 


-4.040 


-4.050 


-4.059 


16 o 


P 


lSl/2 


-28.7634 


-28.7722 


-28.7810 


-28.7897 


-28.7985 




1P3/2 


-13.6182 


-13.6293 


-13.6405 


-13.6516 


-13.6627 




lPl/2 


-11.8432 


-11.8547 


-11.8661 


-11.8775 


-11.8889 


11 


lSl/2 


-32.7665 


-32.7761 


-32.7857 


-32.7952 


-32.8048 




1P3/2 


-17.3721 


-17.3841 


-17.3962 


-17.4082 


-17.4202 




lPl/2 


-15.5826 


-15.5949 


-15.6073 


-15.6196 


-15.6319 




E/nucleon 


-5.941 


-5.951 


-5.961 


-5.971 


-5.981 



the method of moments to the second energy sheet is not 
as simple, because the pole is along the negative imag- 
inary momentum axis. However, since this anti-bound 
state pole is close to the zero energy {E P = —0.066 MeV), 
we have chosen the zero energy point to reduce the 
Schrodinger equation using the method of moments to 
a set of n algebraic equations. It has been demon- 
strated [I?! that this procedure gives a good represen- 
tation of the original potential for the low energy scat- 
tering parameters. As a result we use the effective range 
expansion to determine the position of the anti-bound 
state pole in the momentum or fc-plane, i.e. we write the 
on-shell amplitude in terms of the phase shifts So as 



7r/x k cot Sq — ik ' 

where \i is the reduced mass, and make use of the effective 
range expansion 

fccot^o = -— + i r s k 2 - P s r 3 k A + • ■ • , (28) 
a s 2 

where P s is the shape parameter, to analytically continue 
the amplitude onto the second energy sheet. Since the 
anti-bound state is close to zero energy (A; ~ — 0.04 i), 
we can truncate the effective range expansion to include 
the fc 4 term. To test the accuracy of this procedure, 



we compare the position of the pole on the second en- 
ergy sheet for the Bryan-Scott potential by truncating at 
k 2 and k 4 term with the result Ep — —0.0711531 and 
—0.0711548 MeV respectively. As a result we have cho- 
sen to truncate the effective range expansion to include 
the fc 4 term. 

The use of the trial function in Eq. (|26p has the added 
advantage of allowing us to construct an equivalent rank 
one separable potential, often referred to as the Uni- 
tary Pole Approximation (UPA), that has identically the 
same deuteron wave function as the original OBE poten- 
tial [45[ . After partial wave expansion, this is of the form 

V$ A (k, k') = g e (k) C f . t g t {k') , (29) 

where the form factors are directly related to the radial 
wave function ipe(r) and the strength of the potential is 
adjusted to insure that the matrix element of the UPA 
and original OBE potential are identical at the energy of 
the pole in the amplitude. The same procedure is applied 
to the 4 So channel. 

Having constructed a rank one separable potential 
equivalent to the OBE potential, we can write the Fad- 
deev equations as a set of coupled one dimensional inte- 
gral equations 48]. If one includes the 1 So and 3 Si- 3 Di 
nucleon-nucleon partial waves only, then the number of 
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TABLE VI. Variation in the position of the anti-bound state 
pole on the second energy sheet, the binding energy of the 
deuteron and triton with changes in hadron mass ran- For 
the Bryan-Scott potential the position of the anti-bound 
state pole is Ep = -7.1155E-02 MeV, deuteron binding en- 
ergy Ed = 2.18365 MeV, while the triton binding energy 
E t = 7.9131 MeV in the UPA. 



H 


m H (MeV) 


SE P 


SEp 
5m jf 


SEt 


TV 


138.7 


2.38E-03 


- 0.0201 


-0.0146 


n 


548.7 


-7.40E-05 


0.0019 


0.0034 




550.0 


-1.99E-03 


-0.1026 


-0.3355 




600.0 


-3.09E-03 


0.0486 


0.0790 


p 


763.0 


1.27E-04 


-0.0295 


-0.0517 


OJ 


782.8 


1.34E-02 


0.0923 


0.2776 


N 


938.92 


2.95E-04 


0.0289 


0.0527 



coupled integral equations reduces to five, and these can 
be solved for the binding energy and wave function of the 
triton [H. 

To examine the variation in the binding energy with 
changes in the mass of the mesons and nucleon, we have 
calculated the slope of the binding energy as function 
of the mass at the value of the mass used in the OBE 
potential. In Table I VI I we present this variation in the 
energy of anti-bound state, the deuteron and triton bind- 
ing energies with respect to the variation in the masses 
of the six bosons included in the OBE potential. We 
have also included the variation in the binding energies 
with changes in the nucleon mass rriN- Here, we note 
that the nucleon mass is present, not only in the kinetic 
energy of the two- and three-body equations, but also in 
the definition of the Bryan-Scott OBE potential. For the 
one pion exchange component, the strength of the poten- 
tial is proportional to (g^NN /2M) 2 which is equivalent 
to {fuNN /mir) 2 had BS used a pseudo-vector coupling 
in the Lagrangian. From the Goldberger-Treiman [50| 
relation we have that 



9-xNN 9a 

oc — 

M U 



(30) 



where /„• is the pion decay constant. Although and f v 
are dependent on the quark mass, the ratio to first order 
is not sensitive to variation in quark mass. This suggests 
that the strength of the one pion exchange component of 
the BS should not change with changes in the nucleon 
mass. Since the rj is part of the same SU(3) octet as 
the pion, one could apply the same argument the rj ex- 
change component of the OBE potential. For the scalar 
(<To and <7i) and vector (p and to) meson exchanges, the 
relative strength of the central, the spin-orbit and the 
tensor component depend on the nucleon mass, and to 
that extent, we have maintained the M dependence of 
the OBE potential for the scalar and vector exchanges. 
From Table PVTl we observe that the variation is largest for 



the do an d u), followed by the variation with the tt, <7i, p 
and N masses, with the variation in the energy with the 
rj mass being minimal. 



A. Total variation in binding energies 



From the detailed results given in Table I VII and the 
earlier results for the variation of the meson and nucleon 
masses with quark mass, we can readily deduce the total 
variation of the deuteron and triton binding energies and 
the energy of the anti-bound state, Ep, with changes in 
the quark mass: 



^£ = _ .912^ 

E D m, 



and 



SE t 

~e7 

SE P 
~Ep~ 



-0.980^ 

TO 



-2.839^ 



(31) 



(32) 



(33) 



The details for these calculations are shown in the Ap- 
pendix. 

The variations of the deuteron and triton binding en- 
ergies given in Eqs. (|3"Tj) and (|32|) . respectively, are com- 
pletely compatible with those reported by Flambaum and 
Wiringa [24[. In particular, the coefficiants on the rhs of 
those equations, namely -0.91 for the deuteron and -0.98 
for the triton, are very close to those reported in Ref. [24| 
for the AV14 potential, namely -0.84 and -0.89. 

On the other hand, for the 1 S'o anti-bound state, with 
energy Ep, there is a significant disagreement. The sign 
reported above for SEp/Ep is negative, whereas a pos- 
itive value was reported in Ref. [24|. Since Dmitriev et 
al. [24j presented an apparently general argument relat- 
ing the change in the deuteron binding to that in the 
energy of the anti-bound state, we re-checked every term 
in our calculation carefully. There is no doubt that our 
result is correct for the model used. We note that, from 
Table II of Flambaum and Wiringa (24|, the individual 
pieces of the Argonne potential do not respect the sup- 
posedly general result of Dmitriev et al. and therefore 
it cannot be a model independent result. We note, in 
particular, that the tensor force plays a significant role 
for the deuteron, whereas it is absent for the 1 So chan- 
nel. Clearly, this difference for the 1 S'o anti-bound state 
will lead to significant changes when one computes the 
effect of a change in quark mass on the reaction rate for 
np — > dj. 



VI. CONCLUSIONS 

We have calculated the variation of the binding en- 
ergy of the deuteron, triton and the 1 Sq anti-bound pole 
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position, as well as the binding energy per nucleon for 
a number of light nuclei, with respect to variations in 
the light (average of u and d) quark mass. The results, 
expressed in terms of a parameter Ka, defined by 

~BE{^ = KA — q ' (34) 

are summarised in Table- TVIll In order to determine these 
coefficients, we first calculated the change with quark 
mass of the mesons used in a typical one-boson-exchange 
treatment of the nucleon-nucleon force. Those results 
were summarised in Table Mil For each nucleus we cal- 
culated the rate of change of the binding energy with 
respect to the mass of each meson and the mass of the 
nucleon itself. The values of Ka were obtained by com- 
bining the latter with the results in Table HIT] 

TABLE VII. Coefficients Ka summarising the rate of varia- 
tion of the binding energies and the So anti-bound state pole 
with respect to quark mass - see Eq. (|34p . 



Nucleus 


K A 


D 


-0.912 


T 


-0.979 


E P 


-2.839 


7 Li 


-2.571 


12 C 


-1.438 


16 Q 


-1.082 



For the deuteron our result, Kd = —0.91, is very close 
to that reported by Flambaum and Wiringa 2A] using 
the AV14 potential, namely —0.84. Similarly for the tri- 
ton, our value K t — —0.89 is very close to their value of 
—0.98. The closeness of these results for two rather differ- 
ent treatments of the NN force lends considerable confi- 
dence in their reliability. However, for the position of the 
anti-bound state our calculation differs considerably 
from that of Ref. (24|, taking the opposite sign. This 
suggests that this quantity may be rather more model 
dependent than has been realized hitherto. 

In the case of light nuclei, the binding energies re- 
ported here were calculated in the quark-meson coupling 
(QMC) model, a relativistic mean-field model that takes 
into account the self-consistent response of the internal 
structure of the nucleon to these mean fields. Through 
the self-consistency, the model yields many-body [5l| or 
equivalently density-dependent interactions (52[. Indeed, 
the density dependent Skyrme forces derived from QMC 
have proven remarkably realistic 53] . The values of Ka 
deduced in this way for 7 Li, 12 C and 16 are reported in 
Eqs. (|2~4| - (j2"5j) . It is interesting that the value obtained 
for 7 Li, namely Ki^i = —2.57, is significantly larger than 
that reported in Ref. [Hj], namely -1.03 (AV14) and 
-1.50 (AV18+UIX). These authors did suggest that the 
uncertainty on the value of K could be as large as a factor 



of two and our value is consistent at that level. Clearly, 
this degree of variation calls for more investigation to see 
whether the model dependence can be reduced. 

Our study of these variations of binding energies with 
quark mass is, of course, motivated by the possible ef- 
fects on big bang nucleosynthesis (BBN). Amongst the 
many challenges there, the sizeable discrepancy in the 
abundance of 7 Li with the latest photon-to-baryon ratio 
(post WMAP) is of particular interest. Figure [3] illus- 
trates the 7 Li abundance calculated using the BBN code 
of Kawano [54] , if one allows only the binding energy of 
the deuteron and the energy of the virtual 1 5'o state to 
change with quark mass. The curves correspond to the 
values of Kd and Kp calculated here (solid line) as well 
as the values used by Berengut et al. [55[ (dashed line). 
The substantial difference in slope means that while a 3% 
shift in Sm q /m q would suffice to reproduce the empiri- 
cal abundance using the values of Berengut et al., with 
our values this would require a huge change in quark 
mass. This simple example illustrates the importance of 
a complete study of the BBN problem including all of the 
consequences of a shift of quark mass within the current 
approach, which we leave for future work. Finally, we 



'Li/H 




-0.01 0.00 0.01 0.02 0.03 

<5m,,/m, ( 

FIG. 3. (Color online) Abundance of 7 Li with respect to 
changes in the quark mass in p(n, j)d calculated in the same 
way as [55j (dashed-red line) and using our results for Kd 
and Ke p (continuous- blue line). 

note that while the variation of the light quark masses 
should be most important, it will also be necessary to 
take into account the effect of a corresponding change in 
the strange quark mass, especially now that the strange 
quark sigma commutator seems to be under control [56j . 
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APPENDIX so that: 



From Table I VII we find the variations of the binding 
energies for the deuteron and triton (_Ej with i = D,T), 
and the position of the pole for the anti-bound state 
Ep; according to changes in the mass of the hadrons 
(m H ): 



E, 



1 



SEj 
^— ' 8m h 

H 



dm 



H 



(35) 



We then relate the variation of the mass of each hadron 
to the variation of the quark mass, as given in Eq. (|20l) : 



Smn 
m H 



v H - 



8m a 



(36) 



5m H = (v H ■ m H ) ^ 1 



(37) 



Combining those results we finally obtain the formula 
that gives rise to the results in Eqs. (|3"TT) . (|32l) and (|33p : 



SEi _ 1 ^ SEi 
E, ~ E t . Sm H 



(v H ■ m H ) * 



(38) 
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